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ABSTRACT 

Gamma-ray burst (GRB) observations at very high energies (VHE, E > 
100 GeV) can impose tight constraints on some GRB emission models. Many 
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GRB afterglow models predict a VHE component similar to that seen in blazars 
and plerions, in which the GRB spectral energy distribution has a double-peaked 
shape extending into the VHE regime. VHE emission coincident with delayed X- 
ray flare emission has also been predicted. GRB follow-up observations have had 
high priority in the observing program at the Whipple 10 m Gamma-ray Tele- 
scope and GRBs will continue to be high priority targets as the next generation 
observatory, VERITAS, comes on-line. Upper limits on the VHE emission, at 
late times (> ~4 hours), from seven GRBs observed with the Whipple Telescope 
are reported here. 

Subject headings: gamma rays: bursts — gamma rays: observations 



1. Introduction 



Since their discovery in 1969 (IKlebesadel. Strong fc Olsenl Il973l ). gamma-ray bursts 
(GRBs) have been well studied at many wavelengths. Although various open questions 
remain on their nature, there is almost universal agreement that the basic mechanism is an 
expanding relativistic fireball, that the radiation is beamed, that the prompt emission is due 
to internal shocks and that the afterglow arises from external shocks. It is likely that Lorentz 
factors of a few hundred are involved, with the radiating particles, either electrons or pro- 
tons, being accelerated to very high energies. GRBs are sub-classified into two categories, 
long and short bursts, based on the timescale over which 90% of the prompt gamma-ray 
emission is detected. 



Recently, the Swift GRB Explorer (jGehrels et al.ll2004f) has reveale d that many GRBs 
have associated X-ray flares ( Burrows et al.ll2005l ; iFalcone et al.ll2006a ). These flares have 
been detected between 10 2 s and 10 5 s after the initial prompt emission and have been found 
to have fluences ranging from a small fraction of, up to a value comparable to, that con- 
tained in the prompt GRB emission. This X-ray flare emission has been postulated to 
arise from a number of different scenarios, including late central engine activity where 
the GRB progenitor remains ac t ive for some time after, or re-activates after, the initial 



explosion (IKumar fc 



Proga fc Zhang 



2006 



iran 



200( 



3 Izhang et all bood : iNousek et all I2OO6I : IPerna et all 12005 



King et al.l 120051 ) and refreshed shocks which occur when slower mov- 



ing shells eject ed by the central engine in the prompt phase catch up with the afterglow shock 



at late times (IRees fc Meszarod 1 19981 : ISari fc Meszarod |2000e iGranot. Nakar fc Piran 1 12003 



Guetta et al.ll2006h. For short GRBs, shock heating of a b inary stellar companion has also 
been proposed (IMacFadyen. Ramirez- Ruiz fc Zhang! 120051 ). It is not yet clear whether the 
X-ray flares are the result of prolonged central engine activity, refreshed shocks or some other 
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mechanism (IPanaitescu et al.ll2006l ). A very high energy (VHE; E > 100 GeV) component of 
this X-ray flare emission has also been predicted (jWang. Li fc Meszarosll2006l ). 



Within the standard fireball shock scenario (IRees fc Meszaroa Il992l ; iMeszaros fc Rees 



19931 ; ISari. Piran fc Narayanlll998l ). many models have been proposed which predict emission 
at and above GeV energies during both the prompt an d afterglow phases of the GRB. These 
have been summarized by IZhang &: Meszarod (120041 ) . and references therein, and include 
leptonic models in which gamma rays are produced by electron self-inverse- Compton emission 
from the internal shocks or from the external forward or reverse shocks. Other models 
predict gamma rays from proton synchrotron emission or photomeson cascade emission in 
the external shock or from a combination of proton synchrotron emission and photomeson 
cascade emission from internal shocks. 

Although GRB observations are an important component of the program at many VHE 
observatories, correlated observations at these short wavelengths remain sparse even though 
tantalizing and inherently very important. The sparsity of observations of GRBs at energies 
above lOMeV is dictated not by lack of interest in such phenomena, or the absence of 
theoretical predictions that the emission should occur, but by the experimental difficulties. 

For the observation of photons of energies above 100 GeV, only ground-based telescopes 
are available at present. These ground-based telescopes fall into two broad categories, air 
shower arrays and atmospheric Cherenkov telescopes (of which the majority are Imaging 
Atmospheric Cherenkov Telescopes, or IACTs). The air shower arrays, which have wide 
fields of view making them particularly suitable for GRB searches, are relative ly insensitive. 



There are several reports from these instruments of possible TeV emission. iPadilla et al. 



(Il998h reported possible VHE emission at E > 16 TeV from GRB 920925c. While finding 
no individual burst which is statistically significant, the Tibet-AS" Y Collaboration found a n 
indication of 10 TeV emission in a stacked analysis of 57 bursts ( lAmenomori et al.ll200ll ). 
The Milagro Collaboration reported on the detection of an exces s gamma-ray signa l during 
the prompt phase of GRB 970417a with the Milagrito detector (jAtkins et al.ll2000l ). In all 
of these cases however, the statistical significance of the detection is not high enough to be 
conclusive. In a ddition to searching the Milagro data for VHE counterparts for 25 satellite- 
triggered GRBs (jAtkins et al.ll2005l ). the Milagro Collaboration con ducted a search for VHE 
transients of 40 seconds to 3 hours duration in the northern sky (jAtkins et al.ll2004l ). No 
evidence for VHE emission was found from either of these searches and upper limits on the 
VHE emission from GRBs were derived. Atmospheric Cherenkov telescopes, particularly 
those that utilize the imaging technique, are inherently more flux-sensitive than air shower 



arrays and have better energy resolution but are limited by their small fields of view 
and low duty cycle (~ 7%). In the Burst And Transient Source Explorer (jMeegan et al 



(3-5°) 



1992( 1 
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era (1991-2000), attempts a t GRB monitoring were l imited by slew times and uncertainty 
in the GRB source position (jConnaughton et al.lll997l ). 



Swift, the first of the next generation of gamma-ray satellites which will include AG- 
ILE (Astro-rivelatore Gamma a Immagini LEggero) and the Gamma-ray Large Area Space 
Telescope (GLAST), is beginning to provide arcminute localizations so that IACTs are no 
longer required to scan a large GRB error box in order to achieve full coverage of the possible 
emission region. The work in this paper covers the time period prior to the launch of the 
Swift satellite. 

The minimum detectable fluence with an IACT, such as the Whipple 10m, in a ten sec- 
ond integration is < 10~ 8 erg cm -2 (5 photons of 300 GeV in 5x 10 8 cm 2 collection area). This 
is a factor of > 100 better than GLAST will achieve (3 photons of 10 GeV in 10 4 cm 2 collection 
area). This ignores the large solid angle advantage of a space telescope and the possible steep- 
ening of the observable spectrum because of the inherent emission mechanism and the effect of 



GRB emission in anc 


above the GeV energy rang 


8 (Meszaros, Rees, & Papathanassiou 1994: 


Boettcher & Dermei 


19981: Pi 


lla & Loebl 19981: 


Wang, Dai & Lu 2 


J001; 


Zhang & Meszaros! 


2001; 


Guetta & Granotl2003bl: 


Dermer & Atovan 


2004: Frasdle et al 


2004 


) ; also see Zhang & Meszaros 


(2004 


) and references therein. 



Until AGILE and GLAST are launched, the GRB observations that were made by 
the EGRET experiment on the Compton Gamma Ray Observatory (CGRO) will remain the 
most constraining in the energy range from 30 MeV to 30 GeV. Although EGRET was limited 
by a small collection area and large dead time for GRB detection, it made sufficient detections 
to indicate that there is a prompt component with a hard spectrum that extends at least to 
100 MeV energies. The average spectrum of four bursts detected by EGRET (GRBs 910503, 
9301 31, 940217 & 940301) did not show any evidence for a cutoff up to 10 GeV JDingus 
200ll ). The relative insensitivity of EGRET was s uch that it was not possible to eli minate the 
possibility that all GRBs had hard components (IDingus. Catelli fc Schneidlll998l ). EGRET 
also detected an afterglow component from GRB 940217 that extended to 18 GeV for at least 
1.5 hours after the prompt emission indicating that a high-energy spectral co mponent can 



exten d into the GeV band for a long period of time, at least for some GRBs ( IHurley et al. 



19941 ) . The spectral slope of this component is sufficiently flat that its detection at still 
higher energies may be possible (jMannheim. Hartmann fc Funk! Il996l ). iMeszaros fc Rees 
(119941 ) attribute this emission to the combination of prompt MeV radiation from internal 
shocks with a more prolonged GeV inverse Compton component from external shocks. It is 
also postulate d that this emission could be the result of inverse Compton scattering of X-ray 
flare photons (jWang. Li Meszaros! 120061 ). Although somewhat extreme parameters must 
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be assumed, synchrotron self-Comp ton emission from t he re verse shock is cited as the best 



candidate for this GeV emission by iGranot &: Guettal (120031 ). given the spectral slope that 



was recorded. This requirement of such extreme parameters nat urally explains the lack o f 



GRBs for which such a high energy component has been observed. iGuetta fc Granotl (12003a! ) 
postulate that some GRB explosions occur inside pulsar wind bubbles. In such scenarios, 
afterglow electrons upscatter pulsar wind bubble photons to higher energies during the early 
afterglow thus producing GeV emission such as that observed in GRB 940217. 

The GRB observational data are extraordinarily complex and there is no complete and 
definitive explanation for the diversity of properties observed. It is important to estab- 
lish whether there is, in general, a VHE component of emission present during either the 
prompt or afterglow phase of the GRB. Understanding the nature of such emission will 
provide important information about the physical conditions of the emission region. One 
definitive observation of the prompt or afterglow emission could significantly influence our 
understanding of the processes at work in GRB emission and its aftermath. 

In this paper, the GRBs observed with the Whipple 10 m Gamma-ray Telescope in 
response to HETE-2 and INTEGRAL notifications are described. The search for VHE emis- 
sion is restricted to times on the order of hours after the GRB. In Section [2J the observing 
strategy, telescope configuration, and data analysis methods used in this paper are described. 
The properties of the GRBs observed and their observation with the Whipple Telescope are 
described in Section [3J Finally, in Section HJ the results are summarized and their implica- 
tions discussed in the context of some theoretical models that predict VHE emission from 
GRBs. The sensitivity of future instruments such as VERITAS to GRBs is also discussed. 



2. The Gamma-ray Burst Observations 
2.1. Telescope Configuration 



The observations presented here were made with the 10 m Gamma-ray Telescope at the 
Fred Lawrence Whipple Observatory. Constru cted in 1968, the telescope has been operated 
as an IACT since 1982 (IKildea et al.llin press! ) . In September 2005, the observing program 
at the 10 m was redefined and the instrument was dedicated solely to the monitoring of 
TeV blazars and the search for VHE emission from GRBs. Located on Mount Hopkins 
approximately 40 km south of Tucson in southern Arizona at an altitude of 2300 m, the 
telescope consists of 248 hexagonal mirror facets mounted on a 10 m spherical dish with an 
imaging camer a at its focus. The fron t- aluminized mirrors are mounted using the Davies- 
Cotton design (IDavies Cotton! 119571 ). 



- 7- 



The imaging camera consists of 379 photo-multiplier tubes (PMTs) arranged in a hexag- 
onal pattern. A plate of light-collecting cones is mounted i n front of the PMTs to increase 
their light-collection efficiency. A pattern-sensitive trigger (IBradbury fc Rose! 120021 ) . gener- 
ates a trigger whenever three adjacent PMTs register a signal above a level preset in the 
constant fraction discriminators. The PMT signals for each triggering event are read out 
and digitized using charge-integrating analog to digital converters. In this way, a map of the 
amount of charge in each PMT across the camera is recorded for each event and stored for 
offline analysis. The telescope triggers at a rate of ~ 25 Hz (including background cosmic 
ray triggers) when pointing at high (> 50°) elevation. Although sensitive in the energy range 
from 200 GeV to 10 TeV, the peak response energy of the telescope to a Crab-like spectrum 
during the observations reported upon here was approximately 400 GeV. This is the energy 
at which the telescope is most efficient at detecting gamma rays and is subject to a 20% 
uncertainty. 



2.2. Observing Strategy 

Burst notifications at the Whipple Telescope for th e observations describ ed here were 



received via email from the Global Coordinates Network (IGCN Webpagdl2006l ). When a no- 
tification email arrived, the GRB location and time were extracted and sent to the telescope 
tracking control computer. An audible alarm sounded to alert the observer of the arrival of 
a burst notification. If at sufficient elevation, the observer approved the observations and 
the telescope was commanded to slew immediately to the location of the GRB. The Whipple 
Telescope slews at a speed of 1° s" 1 and therefore can reach any part of the visible sky within 
three minutes. 

Seven different GRB locations were observed with the Whipple 10 m Telescope between 
November 2002 and April 2004. These observations are summarized in Table [TJ At the time 
these data were taken, the point spread function of the Whipple Telescope was approximately 
0.1° which corresponds to the field of view of one PMT. The positional offsets for the GRB 
observations (see Table [2]) were all less than this so a conventional "point source" analysis 
was performed. 



2.3. Data Analysis 



The data were analyzed using the imaging; technique and analys is procedures pioneered 
and developed by the Whipple Collaboration ([Reynolds et al.lll993l ). In this method, each 
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image is first cleaned to exclude the signals from any pixels that are most likely the result 
of noise. The cleaned images are then characterized by calculating and storing the first, 
second, and third moments of the ligh t distribution in each image. The parameters and 



this procedure are described elsewhere ( [Reynolds et al.lll993l ). Since gamma-ray images are 



known to be compact and elliptical in shape, while those generated by cosmic ray showers 
tend to be broader with more fluctuations, cuts can be derived on the above parameters 
which reject approximately 99.7% of the background images while retaining over 50% of 
those generated by gamma-ray showers. These cuts are optimized using data taken on the 
Crab Nebula which is used as the standard candle in the TeV sky. 

Two differe nt modes of observati on are employed at the Whipple Telescope, 11 On - Off 1 
and "Tracking" ( Catanese et al.|[l998l ). The choice of mode depends upon the nature of the 



target. The GRB data presented here were all taken in the Tracking mode. Unlike data taken 
in the On - Off' mode, scans taken in the Tracking mode do not have independent control 
data which can be used to establish the background level of gamma-ray like events during 
the scan. These control data are essential in order to estimate the number of events passing 
all cuts which would have been detected during the scan in the absence of the candidate 
gamma-ray source. In ord er to perform this estimate, a tracking ratio is calculated by 



analyzing "darkfield data" (IHoran et al.ll2002l ). These consist of Off-source data taken in 



the On - Off mode and of observations of objects found not to be sources of gamma rays. 
A large database of these scans is analyzed and in this way, the background level of events 
passing all gamma-ray selection criteria can be characterized as a function of zenith angle. 
Since the GRB data described in this paper were taken at elevations between 50° and 80°, 
a large sample of darkfield data (~ 233 hours) spanning a similar zenith angle range was 
analyzed so that the background during the gamma-ray burst data runs could be estimated. 



3. The Gamma-ray Bursts 

This paper concentrates on the GRB observations made in response to HETE-2 and 
INTEGRAL triggers with the Whipple 10 m Gamma-ray T elescope; observations made in 



response to Swift triggers are the subject of a separate paper (IDowdall et al.llin prep.l ). When 
the GRB data were filtered to remove observations made at large zenith angles, during in- 
ferior weather conditions, and of positions later reported to be the result of false triggers 
or to have large positional errors, the data from observations of seven GRB locations re- 
mained. These GRBs took place between UT dates 021112 and 040422; two have redshifts 
derived from spectral measurements, one has an estimated redshift and four lie at unknown 
distances. Five of the sets of GRB follow-up observations were carried out in response to 
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triggers from the high energy transient explorer 2 (HETE-2; lLamb et al.l (120001 )) while two 
sets of observ ations were triggered by the international gamma-ray astrophysical laboratory 
(INTEGRAL; I Winkler et al.l (119991 )). In the remainder of this section, the properties of each 
of the GRBs observed and the results of these observations are presented. A summary of the 
GRB properties is given in Tabled] while the observations taken at the Whipple Observatory 
are summarized in Table [2j 



3.1. GRB 021112 



This was a long GRB with a durat ion of > 5 s and a peak flux of > 3 x 10 8 erg cm 2 s 



in 



the 8-40 keV band (IRicker et al.ll2002l ) . In the 30 - 400 keV energy band, the burst had a peak 



energ y of 57.15 keV, a duration of 6.39s and a fluence of 2.1 x 10 7 erg cm 2 (IMIT Webpages 



20061 ). The triggering instrument was the French Gamma Telescope (FREGATE) instru- 
ment on HETE-2. The Milagro data taken during the time of this burst were searched for 
GeV/TeV gamma-ray emission. No evidence for prompt emission was found and a prelim- 
inary analysis, assuming a differential photon spectral index of -2.4, gave an upper limit 
on the fluence a t the 99.9% confidence level of J(0.2- 20 TeV) < 2.6xl0 _6 ergcm -2 over a 5 



second interval (IMcEnery et al.ll2002al ). Optical observations with the 0.6-meter Red Buttes 
Observatory Telescope beginning 1.8 hours after the burst did not show any evidence for 
an optical counterpart and placed a limiting magnitude of i? c =21.8 (3sigma) on the optical 
emission: at the time, this was the de epest non-detection of an optical afterglow within 2.6 
hours of a GRB JSchaefer et al.ll2002h . 



Two sets of observations on the location of GRB 021112 were made with the Whipple 
10 m Telescope. The first observations commenced 4.2 hours after the GRB occurred and 
lasted for 110.6 minutes. Observations were also taken for 55.3 minutes on the following 
night, 28.6 hours after the GRB occurred. Upper limits (99.7% c.l.) of 0.20 CratQ and 0.30 
Crab (E > 400 GeV), respectively, were derived for these observations assuming a Crab-like 
spectrum (spectral index of -2.49). 



1 Since the Crab is the standard candle in the VHE regime, it is customary to quote upper limits as a 
fraction of the Crab flux at the same energy. 
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3.2. GRB 021204 



Little information is available in the literature on this HETE-2 burst. The GRB location 



Schmidt at Kiso Observatory (jUrata et al. 



was o bserved with a number of o ptical telescopes (the RIKEN 0.2m (iTorii, Yamaoka fe Kato 
20021 ) . the 32 inch Tenagra I I (INvsewander. Reichart fe Schwarta 120021 ). and the 1.05 m 



20021 )) but no optical transient was found to a 



limiting magnitude of i?=16.5, 2.1 ho urs after the burst (jTorii. Yamaoka fe Katoll2002l ). and 
to #=18.8, 6.2 hours after the burst JUrata et alJbooj ). 



Whipple observations of this burst location commenced 16.9 hours after the GRB oc- 
curred and lasted for 55.3 minutes. An upper limit (99.7% c.l.) of 0.33 Crab was derived for 
the VHE emission above 400 GeV during these observations. 



3.3. GRB 021211 



This long, bright burst was detected by all three instruments on HETE-2. It had a 
duration > 5.7s in the 8-40keV band with a fluence of 

Y-7 



10 6 erg cm 2 during that in terval 



(I Crew 2002 ). The peak flux was >8xl0 7 erg cm 2 s 1 (i.e. > 20 Crab flux) in 5 ms (I Crew 
2002h . This burst had a peak energy of 45.56 ke V, a duration of 2.80s , and a fluence o f 



2.4xl0" 6 ergcm- 2 in the 30-400keV energy band flMIT Webpagesl 120061 ). iFox et all (12003 ) 



reported on the early optical, near-infrared, and radio observations of this burst. They 
identified a break in the optical light curve of the burst at t=0.1-0.2hr, which was inter- 
preted as the signature of a reverse shock. The light curve comprised two distinct phases. 
The initial steeply-declining flash was followed by emission declining as a typical afterglow 
with a power-law index close to 1. KAIT observations of the afterglow also detected the 



2002: Lamb et al. 2002; Garnavich et al. 2002a 



steeply declining light curve and evidence for an early break (|Li et al.l 12003). The optical 
trans i ent was detected at many observatories ( 



'ark. Williams fe Barthelmy 



2002: Li et al. 



b|). The optical transient faded from an R 



band magnitude of 18. 3, 20.7 minutes af ter the burst, to an R- band magnitude of 21.1, 5.7 
hours after the burst (IFox et al.l 120031 ) . IVreeswijk et al.l (120021 ) derived a redshift of 1.006 
for this burst based on spectroscopic observations carried out with the European Southern 
Observatory's Very Large Telescope (VLT) at Paranal, Chile. Milagro searched for emis- 
sion at GeV/TeV energies over the burst duration reported by the HETE-2 wide field X-ray 
monitor. They did not find any evidence for prompt emission and a preliminary analysis, 
assuming a differential photon spectral index of -2.4, gave an upper limit on the fluence 
at the 99.9% confidenc e level of J(0.2- 20 TeV)<3.8x 10~ 6 erg cm -2 over a 6 second interval 
jMcEnerv et al.ll2002bh . 
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Whipple observations on this GRB location were initiated 20.7 hours after the GRB 
and lasted for 82.8 minutes. An upper limit (99.7% c.l.) on the VHE emission of 0.33 Crab 
(E > 400 GeV) was derived from these observations. 



3.4. GRB 030329 



This GRB is one of the brightest bursts on record. It triggered the FREGATE instru- 
ment on HETE-2 in the 6- 120 keV energy band. It had a duration 22.76 seco nds, a ffuence of 



l.lx 10 4 erg cm 2 and a peak energy of 67.86 keV in the 30-400keV band (IMIT Webpages 



20061 ). The peak flux over 1.2 seconds was 7x10 6 er g cm 
Crab flux in that energy band (IVanderspek et al.l 120031 ) . 



: s 1 which is > 100 times the 



The optica l transient was identified by lPeterson fc Pried (120031 ) . Due to its slow decay 
(jUemura 1120031 ) and brightness (i?~13), extensive photometric observations were possible, 
making th is one of the best-o bserved GRB afterglows to date. Early observations with 
the VLT (IGreiner et al.l 120031 ) revealed evidence for narrow emission lines from the host 
galaxy indicating that this GRB occurred at a low redshift of z= 0.1687. Observations 
of the afterglow continued for many nights as it remained bright with a slow but uneven 
rate of decline and exhibited some episodes of increasing brightness. These observations 
are well-documented in the GCN archives. Spectral measurements made on 6 April 2003 



by IStanek et al.l (j2003al ) showed the development of broad peaks in flux, characteristic of 
a supernova. Over the next few nights , the afterglow emis sion faded and the features of 
the supernova became more prominent ( IStanek et al.l l2003bl ) . These observations provided 
the first direct spectroscopic evidence that at least a subset of GRBs is associated with 
super novae. 

The afterglow was detected at many other wavelengths. Radio observations with the 
VLA det ected a 3.5 mJy source at 8.46 G Hz. This is the brightest radio afterglow detected 
to date (IBerger. Soderberg &: Fraill 120031 ) . The aft erglow was also b right at sub millimeter 
(jSmith et al.ll2003l ) and near infrared wavelengths (ILamb et al.ll2003l ). The X-ray afterglow 
was detect ed by RXTE during a 27 -minute observation that began 4 hours 51 minutes after 



the burst ([Marshall fc Swankll2003l ). The flux was ~1.4x 10" 10 erg cm" 2 s' 1 in the 2-10keV 



band (-0.007% of the Crab). 

Whipple observations of the location of GRB 030329 commenced 64.6 hours after the 
prompt emission. In total, 241.4 minutes of observation were taken spanning five nights. 
The upper limits (99.7% c.l.) from each night of observation are listed in Tableland are 
displayed on the same temporal scale as the optical light curve of the GRB afterglow in 



Figured] When these data were combined, an upper limit (99.7% c.l.) for the VHE emission 
above 400 GeV of 0.17 Crab was derived. 



3.5. GRB 030501 



This burst was initially detected by the imager on board t he INTEGRAL satellit e 
(IBIS/ISGRI) and was found to have a duration of ~40 seconds (jMereghetti et al.l 120031 ) . 
The burst was also detecte d by the Ulysses spacecraft and the spectrometer instrument 
(SPI-ACS) on INTEGRAL ( iHurley et al.ll2003l ). Triangulation between these two detections 
allowed a position annulus to be computed for this GRB. As observed by Ulysses, it had a 
duration of 

with a peak flux of 4.9x10 



75 seconds and had a 25- 100 keV fluence of approximately l.lx 10 6 erg cm 2 



" 7 erg cm 2 s 1 over 0.25 seconds. Follow-up optical observa- 



i?=20.0, 16.5 hours after the burst (jOfek et al 



i?=18.0, 0.3- 17 minutes after the b urst (Boer fc Klotj 20031 ) and to a limiting magnitude of 



tions with several telescopes did not find ev i dence for an optical transient (jOfek et al.l 12003 
Rumyantsev. Pavlenko &: Pozanenkdl2003at iBoer fc Klota 120031 ) to a limiting magnitude of 



2003 



Whipple observations of this burst location commenced 6.6 hours after its occurrence 
and continued for 83.1 minutes. An upper limit (99.7% c.l.) on the VHE emission (E > 
400 GeV) during these observations of 0.27 Crab was derived. 



3.6. GRB 031026 



This burst was located by the FREGATE instrument on HETE-2. It had a dura- 
tion of 114.2 second s with a fluence of 2.3x 10~ 6 ergcm -2 in the 25-100keV energy band 
(IRicker et al.ll2003al ) while in t he 30-400keV energy band it had a duration of 31.97 s and 



a fluence of 2.8x10 6 erg cm 2 ( 1MIT Webpaged 120061 ) . Follow-up optical observations were 



carried out with a numbe r of instruments including the 1.05 m Schmidt at the Kiso O bser- 
vatory (jUrata et al.l 120031 ). the 32 inch Tenagr a II Telescope (jNvsewander et al.ll2003al) . and 
the 1.0m Telescope at the Lulin Observatory (IHuang et al.ll2003al ). but no optical transient 
was fo und to a limiting m agnitude of i?=20.9 for observations taken 6-12 hours after the 



burst (IHuang et al. 



2003a) and to J r =20.4 from observations taken 3.9 and 25.7 hours af- 



ter the burst (jNvsewander et al.ll2003bl ). The 30m IRAM Telescope was used to search the 
field aroun d the GRB location but did not detect any source with a 250 GHz flux density 
> 16mJy (IBertoldi et al.l 120031 ). A spectral analysis of the prompt X-ray and gamma-ray 
emission from this burst revealed it to have a very hard spectrum which is unusual for such 
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a long and relatively faint burst (IRicker et al.ll2003bl ). It was noted that the counts ratio of 
>1.8 between the 7-30 keV and 7 - 80 keV FREGATE energy bands was one of the most ex- 
treme measured (IRicker et al.ll2003al ). A "new pseudo-redshif t" of 6.67 ±2.9 was computed 
for this burst using the prescription of iPelangeon et al.l (120061 ) . 



Whipple observations of this burst location were initiated 3.7 minutes after receiving 
the GRB notification. The burst notification however, was not received until more than 3 
hours after the prompt GRB emission. Although Whipple observations commenced 3.3 hours 
after the prompt emission, the first data run is not included here due to inferior weather 
conditions. The data presented here commenced 3.7 hours after the GRB and continued for 
82.7 minutes. An upper limit (99.7% c.l.) for the VHE emission (£>400GeV) of 0.41 Crab 
was derived. 



3.7. GRB 040422 



This burst was detected by the imaging instrument (IBIS/ISGRI) on the INTEGRAL 
satellite in the 15-200 keV energy band. It had a duration of 8 seconds, a peak flux 
between 20 and 20 keV of 2.7 phot ons cm" 2 s -1 and a fluence (Is integration time) of 



2.5x10 erg cm (IGotz et al.l 120031 ) . Foll ow-up observations were carried out by many 



20031: 



groups but no optical transient was detected (Male sani et al. 2003; ISonoda. Maeno fc Yamauchil 
" llRvkofjl2003l : Irluang et al.1l2003bl : IPiccioni et al.ll2003l : lRumvantsev & Pozanenkoll2003bl : 



Qiu fc Hull2003l ). The ROTSE-IIIb Telescope at McDonald Observatory began taking unfil- 
tered optical data 22.1s after the GRB. Using the first 110 s of data, a li miting magn itude 
of ~ 17.5 was placed on the R-band emission from the GRB at this time (IRykoffl 120031 ) . 



Whipple observations of this burst commenced 4.0 hours after the prompt emission and 
continued for 27.6 minutes. An upper limit (99.7% c.l.) on the VHE emission (E > 400 GeV) 
of 0.62 Crab was derived. 



4. Results and Discussion 

Upper limits on the VHE emission from the locations of seven GRBs have been derived 
over different timescales. For each GRB, a number (1-10) of follow-up 28-minute duration 
observations were taken with the Whipple 10 m Telescope. These GRB data were grouped 
by UT day and were combined to give one upper limit for each day of observation. The 
limits range from 20% to 62% of the Crab flux above 400 GeV and are presented in Table |5J 
In addition to calculating upper limits on the GRB emission for each day, upper limits 
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were calculated for each of the 28-minute scans. These are plotted for each of the GRBs in 
Figure [21 

The usefulness of the upper limits presented here is limited by the fact that five of the 
GRBs occurred at unmeasured redshifts thus making it impossible to infer the effects of the 
infrared background light on those observations. In addition to this, the earliest observation 
was not made at Whipple until 3.68 hours after the prompt GRB emission. Although the 
Whipple 10 m Telescope is capable of beginning GRB observations less than 2 minutes after 
receiving notification, a number of factors, including notifications arriving during daylight 
and delays in the distribution of the GRB locations, delayed the commencement of the GRB 
observations presented here. Although data-taking for GRB 031026 began 3.7 minutes after 
the GRB notification was received, this notification was not distributed by the GCN until 
3.3 hours after the GRB had occurred. Thus, the observations presented here cannot be 
used to place constraints on the VHE component of the initial prompt GRB emission and 
pertain only to the afterglow emission and delayed prompt emission from GRBs. 

One of the main obstacles for VHE observations of GRBs is the distance scale. Pair 
production interactions of gamma rays with the infrared photons of the extragalactic back- 
ground light attenuate the gamma-ray signal thus limiting the distance over which VHE 
gamma rays can pro pagate. Recently howe ver, the H.E.S.S. telescopes have detected the 
blazar PG 1553+113 (lAharonian et al.ll2006l ). The redshift of this object is not known but 
there are strong indications that it lies at z>0.25, possibly as far away as z=0.74. This 
could represent a large increase in distance to the most distant detected TeV source, re- 
vealing more of the universe to be visible to TeV astronomers than was previously thought. 
Although GRBs lie at cosmological distances, many have been detected at redshifts accessible 
to VHE observers. Of the GRBs studied 
while the redshift of one wa s estimated by 



rere, only 2 had spectr oscopic redshifts measured 



Pel angeon et al.l (120061 ) using an improved version 



of the redshift estimator of lAtteia et al.l (120031 ). Since all of the GRBs discussed here were 
long bursts, it is likely that their redshifts are of order 1. Due to the unknown redshifts of 
most of the bursts and the uncertainty in the density of the extragalactic background light, 
the effects of the absorption of VHE gamma rays by the infrared background light have not 
been included here. 



Granot. Nakar fc Piran I (120031 ) analyzed the late time light curve of GRB 030329 and 
find that the large variability observed at several times (t=1.3-~1.7 days, ~ 2.4 -2.8 days, 
~ 3.1 -3.5 days and at ~ 4.9 -5.7 days) after the burst is most likely the result of refreshed 
shocks. These time intervals have been highlighted in the top panel of Figure [T] and it can 
be seen that some of the observations taken at Whipple occurred during these times thus 
imposing upper limits on the VHE emission during these refreshed shocks. Since GRB 030329 
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occurred at a low redshift (z =0.1685), it is possible that the effects of infrared absorption 
on any VHE emission component may not have been significant enough to absorb all VHE 
photons over the energy range to which Whipple is sensitive. 

Figure [3] shows these scan-by-scan upper limits as a function of time since the prompt 
GRB emission. Also plo tted are the predicted fluxe s at various times aft er the GRB by 
Zhang fc Meszarosl (boOlh and lPe'er fc WaxmarJ tooi ) at ~ 400 GeV, and bv lGuetta k Granot 
(j2003al ) at 250 GeV. Although the peak response energy of the Whipple Telescope at the 
time of these observations was 400 GeV, it still had sensitivity, albeit somewhat reduced, at 
250 GeV. 



Razzaque. Meszaros fc Zhang! (120041 ) predict a delayed GeV component in the GRB 
afterglow phase from the inverse-Compton up-scattering on external shock electrons. The 
duration of such a comp onent is predicted to be up to a few hours, softening with time. 
Zhang fc Meszarosl (120011 ) investigated the different radiation mechanisms in GRB afterglows 
and identified parameter-space regimes in which different spectral components dominate. 
They found that the inverse-Compton GeV photon component is likely to be significantly 
more important than a possible proton synchrotron or ele c tron s ynchrotron component at 
these high energies. The predictions of IZhang fc Meszarosl (120011 ) for VHE emission at dif- 
ferent times after a typical "Regime II" burst are shown by squares on Figure [3j Although 
the observations presented here do not constrain these predictions, the sensitivity is close to 
that required to detect the emission predicted. 

A recent analysis of archival data from the EGRET calorimeter has found a multi-MeV 
spectral compon ent in the prompt pha se of GRB 941017, that is distinct from the lower en- 



ergy component (j Gonzalez et al. 



20031 ). This high energy component appeared between 10 s 



and 20 s after the start of the GRB and had a roughly constant flux with a relatively hard 
spectral slope for ~ 200 s. This observation is difficult to explain within the standard syn- 
chrotron model, thus indicating the existence of new phenomena. iGranot fc Guettal (120031 ) 
investigated possible scenarios for this high energy spectral component and found that most 
models fail. They concluded that the best candidate for the emission mechanism is syn- 
chrotron self-Compton emission from the reverse shock and predicted that a bright optical 
transien t, similar to that observe d in GRB 990123, should accompany this high energy com- 
ponent. iPe'er fc Waxmanl (120041 ) explain this high energy tail as emission from the forward 
shock electrons in the early afterglow phase. These electrons inverse-Compton scatter the 
optical photons that are emitted by the reverse shock electrons resulting in powerful VHE 
emission for 100 s to 200 s after the burst as indicated by the lines on Figure El Although the 
observations presented here did not commence early enough after the prompt GRB emission 
to constrain such models, the sensitivity of the Whipple Telescope is such that the VHE 
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emission predicted by these models would be easily detectable for low redshift bursts. 



The prediction of iGuetta &: Granotl (j2003al ) for VHE emission 5 x 10 3 s after the burst 
from the combination of external Compton emission (the relativistic electrons behind the 
afterglow shock upscatter the plerion radiation) and synchrotron self-Compton emission (the 
electrons accelerated in the afterglow emit synchrotron emission and then upscatter this 
emission to the VHE regime) is indicated by a star on Figure [31 The emission is predicted 
to have a cutoff at ~ 250 GeV due to pair production of the high energy photons with the 
radiation field of the pulsar wind bubble. For afterglows with an external density similar to 
that of the inter-stellar medium, photons of up to 1 TeV are possible. It can be seen that, 



althoug h the upper limits presented here are below the predicted flux from lGuetta fc Granot 



( I2003al ). the observations at Whipple took place after this emission was predicted to have 
occurred. Had data taking at Whipple commenced earlier, the emission predicted by these 
authors should have been detectable for nearby GRBs. 



Razzaque. Meszaros fc Zhang! (120041 ) investigated the interactions of GeV and higher 



energy photons in GRB fireballs and their surroundings for the prompt phase of the GRB. 
They predict that high energy photons escaping from the fireball will interact with infrared 
and microwave background photons to produce delayed secondary photons in the GeV -TeV 
range. Although observations of the prompt phase of GRBs are difficult with IACTs since 
they are pointed instruments with small fields of view which must therefore be slewed to 
respond to a burst notification, observations in time to detect the delayed emission are 
possible. 

There are many emission models which predict significant VHE emission during the 
afterglow phase of a GRB either related to the afterglow emission itsel f or as a VHE compo- 
nent of the X-ray flares that have been observed in many Swift bursts. lO'Brien et al.l (120061 ) 
analyzed 40 Swift bursts which had narrow-field instrument data within 10 minutes of the 
trigger and found that ~ 50% had late (t > T90) X-ray flares. If the bulk of the radiation 
comes via synchrotron radiation as is usually supposed, then by analogy with other sys- 
tems with similar properties (supernova remnants, active galactic nuclei jets), it is natural 
to suppose that there must also be an inverse Compton component b y which photons are 
boosted into the GeV -TeV energy range. This process is described by lPilla fc Loebl (119981 ) 
who discuss the relationship between the energy at which the high energy cutoff occurs, the 
bulk Lorentz factor and the size of the emission region. A high energy emission component 
due to inver s e Com pton emission has also been considered in detail for GRB afterglows by 
Sari fc Esinl (120011 ) ; the predicted flux at GeV - TeV energies is comparable to that near the 
peak of the radiation in the afterglow synchro t ron sp ectrum. Only direct observations can 
confirm whether this is so. IGuetta &: Granotl (l2003bl ) predict that the ~300 GeV photons 
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from the prompt GRB phase will interact with background IR photons, making delayed high 
energy emission undetectable unless the intergalactic magnetic fields are extremely small. 



The Swift GRB Explorer has shown that 



(Burrows et al. 


2005) 


( Connauehtonl 


2002) 



50% of GRBs have one or more X-ray 
i (~ 28 hours) after the prompt emission 



(120061 ) have predicted VHE emission coincident in time with the X-ray flare photons. In this 
model, if the X-ray flares are caused by late central engine activity, the VHE photons are 
produced from inverse Compton scattering of the X-ray flare photons from forward shock 
electrons. If the X-ray flares originate in the external shock, VHE photons can be produced 
from synchrotron self-Compton emission of the X-ray flare photons with the electrons which 
produced them. Should VHE emission be detected from a GRB coincident with X-ray flares, 
the time profile of the VHE emission could be used to distinguish between these two origins 
of the X-ray flares. 

No evidence for delayed VHE gamma-ray emission was seen from any of the GRB 
locations observed here and upper limits have been placed on the VHE emission at various 
times after the prompt GRB emission. Although there are no reports of the detection of 
X-ray flares or delayed X-ray emission from any of these GRBs, it is likely that such emission 
was present in at least some of them given the frequency with which it has been detected 
in GRBs observed by Swift. Indeed, the light curve of GRB 030329 shows large variability 
amplitude a few days after the burst and, as shown in Figure [TJ, Whipple observations were 
taken during these episodes. Apart from this, a measured redshift is only available for one of 
the other bursts observed here and it is possible that the remaining five occurred at distances 
too large to be detectable in the VHE regime. 



Soderberg et al.l (120041 ) reported on an unusual GRB (GRB 031203) that was much less 
energetic than average. Its similarity, in terms of brightness, to an earlier GRB (GRB 980425) 
suggests that the nearest and most common GRB eve nts have not been detected up until now 
because GRB detectors were not sensitive enough (jSazanov. Lutovinov &: Sunyaevi 12004 ) . 
Most GRBs that have been studied up until now lie at cosmological distances. They generate 
a highly collimated beam of gamma rays ensuring that they are powerful enough to be 
detectable at large distances. Both of the less powerful GRBs detected to date occurred 
at conside rably lower reds h ifts; G RB 980425 at 2=0.0085 and GRB 031203 at z= 0.1055. 
Although ISoderberg et al.l (120041 ) conclude that up until now, GRB detectors have only 
detected the brightest GRBs and that the nearest and most com mon GRB events have been 
missed because they are less highly collimated and energetic, iRamirez-Ruiz et al.l (120051 ) 
argue that the observations of GRB 031203 can indeed be the result of off-axis viewing of a 
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typical, powerful GRB with a jet. Should future observations prove there to be a closer, less 
powerful population of GRBs, these would be prime targets for IACTs. 

In the past year, the Whipple Observatory 10m Telescope has been used to carry out 
follow-up observations on a number of GRBs detected by the Sw ift GRB Explorer. Th e 



analysis of these observations will be the subject of a separate paper (IDowdall et al. 



m prep 



The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is currently 
under construction at the Fred Lawrence Whipple Observatory in Southern Arizona. Two 
of the four telescopes are fully operational and it is anticipated that the four-telescope array 
will be operational by the end of 2006. GRB observations will receive high priority and, 
when a GRB notification is received, their rapid follow-up will take precedence over all other 
observations. The VERITAS Telescopes can slew at 1° s" 1 thus enabling them to reach any 
part of the visible sky in less than 3 minutes. When an acceptable (i.e. at high enough 
elevation) GRB notification is received during observing at VERITAS, an alarm sounds to 
alert the observer that a GRB position has arrived. Upon receiving authorization from the 
observer, the telescope slews immediately to the position and data-taking begins. Given that 
the maximum time to slew to a GRB is 3 minutes, and that Swift notifications can arrive 
within 30 s of the GRB, it is possible that VERITAS observations could begin as rapidly as 
2-4 minutes after the GRB, depending on its location with respect to the previous VERITAS 
target. 

As has been shown above, the Whipple 10 m Telescope is sensitive enough to detect 
the GRB afterglow emission predicted by many authors. With its improved background 
rejection and greater energy range, VERITAS will be significantly more sensitive for GRB 
observations than the Whipple 10m Telescope. The VERITAS sensitivity for observations 
of different durations is shown in Figure 0} Based on the assumed rate of Swift detections 
(100 year -1 ), the fraction of sky available to VERITAS, the duty-cycle at its site and the 
sun avoidance pointing of Swift which maximizes its overlap with nighttime observations, it 
is anticipated that ~ 10 Swift GRBs will be observable each year with VERITAS. 
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Table 1 


The propert 


ies of the 


gamma-ray bursts described 


in this work. 




LHsrnvprv 


r Tri crcrpr 




Fluence a 


T90 b 


i—iLlXjL £L V v ) < V \-\-\. \- 


PRT) 


satellite 


Number 


z 


(erg cm -2 ) 


(b) 


[Ke\ ) 


021112 


HETE-2 


2448 




2.1xl0" 7 


6.39 


30-400 


021204 


HETE-2 


2486 










021211 


HETE-2 


2493 


1.006 c 


2.4xl0- 6 


2.80 


30-400 


030329 


HETE-2 


2652 


0.17 d 


l.lxlO" 4 


22.76 


30-400 


030501 


INTEGRAL 


596 




l.lxl0~ 6 


~75 c 


25-100 


031026 


HETE-2 


2882 


6.67 f 


2.8xl0~ 6 


31.97 


30-400 


040422 


INTEGRAL 


1758 




g 


8 h 





a The fluence, where available, is quoted for the energy range given in column 7 
over the duration listed in column 6. For most HETE-2 bursts, this was found 
at: 



http: / / space.mit.edu/HETE/Bursts/Data 



b Except for where a footnote is referenced, the durations in this column are 
T90, the time interval during which 90% of the GRB photons were detected in 
the 30-400 keV energy band. 



Vreeswijk et all (120021 ) . 



Greiner et all (120031 ). 



e The fluence and duration given in the table are from burst observations with 
the Ulysses satellite and the SPI-ACS instrument on INTEGRAL. The event 
was quite weak so there is a factor of 2 uncertainty in the numbers quoted 
iHurlev et alil2003h . Observations with the IBIS/ISGRI inst rument on INTE- 
GRA L alone gave a duration of ~40 seconds for the burst (lMereghetti et al 
2003h . 



f This redshift was determined using the redshift estimator described in 
Pelangeon et al.l (120061 ). 



g The fluence was not quoted for this burst over its 8 secon d duration. It ha d 
a fluence of 2.5 x 10" 



" 7 erg cm 2 when integrated over 1 second ()Gotz et al.ll2003l ) 
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h It was not stated by lGotz et al.l (120031 ) whether or not this duration is T90. 
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Table 2. The VHE GRB observations. 



KjSXD 


rp rp a 

1 GRB - -1- OBS 


Exposure 
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rosition Unset 
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rp rp c 

J- GRB - i-UL 

(hrj 


r lux 

l^uraD ) 


021112 
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110.56 


0.013 


5.1 


< 0.200 




28.63 


55.28 


0.013 


29.0 


< 0.303 


021204 


16.91 


55.34 


0.009 


17.4 


< 0.331 


021211 


20.69 


82.79 


0.058 


21.9 


< 0.325 


030329 


64.55 


65.21 


0.060 


66.2 


< 0.360 




112.58 


83.17 


0.022 


113.8 


< 0.279 




136.23 


37.55 


0.022 


137.0 


< 0.323 




162.14 


27.74 


0.022 


162.4 


< 0.519 




186.16 


27.73 


0.022 


186.4 


< 0.399 


030501 


6.58 


83.10 


0.001 


7.3 


< 0.265 


031026 


3.68 


82.70 


0.007 


4.9 


< 0.406 


040422 


3.99 


27.63 


0.062 


4.2 


< 0.620 



a The time in hours between the start of the GRB and the beginning of 
observations with the Whipple 10m Telescope. 

b The angular separation between the position at which these data were 
taken and the refined location of the GRB. 

c The length of time after the GRB for which the upper limits (ULs) are 
quoted. Since all data are combined to compute the upper limit, the mean time 
of the observations is quoted as the time to which the upper limit pertains. 

d This is the flux upper limit in units of equivalent Crab flux above the peak 
response energy of ~ 400 GeV. Above this energy, the integral Crab flux is 
8.412 ± 1.840 x 10" 11 cm" 2 s _1 . 
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Fig. 1. — Top: The flux upper limits above 400 GeV (99.7% c.l.) on the VHE emission 
from GRB 030329. The time p eriods during which the four bumps in the lightcurve occur 
(IGranot. Nakar fc Piran 1 120031 ) ar e shown as shaded rectangles. Bottom: The optical light 
curve of GRB 030329 taken from iLipkin et al.l (120041 ) . The time since the GRB is shown 
with the same scale on the x-axis of both plots. 
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Fig. 2. — For each GRB location observed, flux upper limits in units of 10 -11 erg cm -2 s _1 
were calculated for each 28-minute scan taken. These are plotted here as a function of the 
time since the GRB prompt emission for each GRB. Only one 28-minute observation was 
made on GRB 040422 so the plot for this GRB is not shown. 
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Fig. 3. — The flux upper limits above 400 GeV for all of the GRBs observed (blue triangles). 
The limits are plotted as a functio n of time since the GRB p rompt emission. The approximate 
flux level at 400 GeV predicted by lPe'er fc Waxmanl (120041 ) is indicated by the red solid lines 
along with the time interval durin g which it is predicted to o ccur; magenta squares show the 
emission at 400 GeV predicted bv lZhang fc Meszaros ( 20011) at various times after the GRB 
prompt emission; the prediction of iGuetta fc Granotl (]2003al ) for VHE emission at 250 GeV 
5 x 10 3 s after the burst from the combination of external Compton and synchrotron self- 
Compton emission is shown by the black star. 
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Fig. 4. — The sensitivity of the VERITAS array for exposures of 50 hours, 5 hours, 0.5 hours, 
and 0.05 hours (i.e., 3 minutes). 



